Global Analysis of DNA Methylation Variation in Adipose Tissue from Twins Reveals Links to Disease-Associated Variants in Distal Regulatory Elements

Epigenetic modifications such as DNA methylation play a key role in gene regulation and disease susceptibility. However, little is known about the genome-wide frequency, localization, and function of methylation variation and how it is regulated by genetic and environmental factors. We utilized the Multiple Tissue Human Expression Resource (MuTHER) and generated Illumina 450K adipose methylome data from 648 twins. We found that individual CpGs had low variance and that variability was suppressed in promoters. We noted that DNA methylation variation was highly heritable (h²_median = 0.34) and that shared environmental effects correlated with metabolic phenotype-associated CpGs. Analysis of methylation quantitative-trait loci (metQTL) revealed that 28% of CpGs were associated with nearby SNPs, and when overlapping them with adipose expression quantitative-trait loci (eQTL) from the same individuals, we found that 6% of the loci played a role in regulating both gene expression and DNA methylation. These associations were bidirectional, but there were pronounced negative associations for promoter CpGs. Integration of metQTL with adipose reference epigenomes and disease associations revealed significant enrichment of metQTL overlapping metabolic-trait or disease loci in enhancers (the strongest effects were for high-density lipoprotein cholesterol and body mass index [BMI]). We followed up with the BMI SNP rs713586, a cg01884057 metQTL that overlaps an enhancer upstream of ADCY3, and used bisulphite sequencing to refine this region. Our results showed widespread population invariability yet sequence dependence on adipose DNA methylation but that incorporating maps of regulatory elements aid in linking CpG variation to gene regulation and disease risk in a tissue-dependent manner.     

Ethanol analysis from biological samples by dual rail robotic autosampler

Detection, identification, and quantitation of ethanol and other low molecular weight volatile compounds in liquid matrices by headspace gas chromatography-flame ionization detection (HS-GC-FID) and headspace gas chromatography-mass spectrometry (HS-GC-MS) are becoming commonly used practices in forensic laboratories. Although it is one of the most frequently utilized procedures, sample preparation is usually done manually. Implementing the use of a dual-rail, programmable autosampler can minimize many of the manual steps in sample preparation. The autosampler is configured so that one rail is used for sample preparation and the other rail is used as a traditional autosampler for sample introduction into the gas chromatograph inlet. The sample preparation rail draws up and sequentially adds a saturated sodium chloride solution and internal standard (0.08%, w/v acetonitrile) to a headspace vial containing a biological sample, a calibrator, or a control. Then, the analytical rail moves the sample to the agitator for incubation, followed by sampling of the headspace for analysis. Using DB-624 capillary columns, the method was validated on a GC-FID and confirmed with a GC-MS. The analytes (ethanol, acetonitrile) and possible interferences (acetaldehyde, methanol, pentane, diethyl ether, acetone, isopropanol, methylene chloride, n-propanol, and isovaleraldehyde) were baseline resolved for both the GC-FID and GC-MS methods. This method demonstrated acceptable linearity from 0 to 1500 mg/dL. The lower limit of quantitation (LOQ) was determined to be 17 mg/dL and the limit of detection was 5 mg/dL.     

Identifying the genetic basis of ecologically and biotechnologically useful functions of the bacterium Burkholderia vietnamiensis

Signature-tagged mutagenesis (STM) was used to identify genetic determinants of fitness associated with two key ecological processes mediated by bacteria. Burkholderia vietnamiensis strain G4 was used as a model bacterium to investigate: phenol degradation as a model of bioremediation, and pea rhizosphere colonization as a prerequisite to biological control and phytoremediation. A total of 1900 mutants were screened and 196 putative fitness mutants identified; the genetic basis of 137 of these mutations was determined by correlation to the G4 genome. The phenol-STM screen was more successful at identifying phenol degradation mutations (83 mutants; 4.4% hit rate) than a conventional agar-based phenol screen (49 mutants, 5319 screened, 0.92% hit rate). The combination of both screens completely defined the components of the TOM pathway in strain G4 and also identified novel accessory genes not previously implicated in phenol utilization. The rhizosphere-STM screen identified 113 mutants (5.9% hit rate); 107 had reduced tag signals indicative of poor rhizosphere colonization (Rhiz-), while six mutants produced high hybridization signals suggesting increased rhizosphere competence (Rhiz+). Competition assays confirmed that 69% of Rhiz- mutants tested (24/35) were severely compromised in their rhizosphere fitness. Seventy Rhiz- mutations mapped to genes with the following putative functions: amino acid biosynthesis (25; 36%), general metabolism (18; 26%), hypothetical (9; 13%), regulatory genes (4; 5.7%), transport and stress (2 each; 2.8% respectively). One of the most interesting discoveries mediated by the rhizosphere-STM screen was the identification of three Rhiz+ mutants inactivated within a single virulence-associated autotransporter adhesin gene; this mutation consistently produced a hyper-colonization phenotype suggesting a highly novel role for this surface adhesin during plant interactions. Our study has shown that STM can be successfully applied to ecologically important microbial interactions, defining the underlying genetic systems important for biotechnological fitness of environmental bacteria such those from the Burkholderia cepacia complex.     

Development of a recA gene-based identification approach for the entire Burkholderia genus

Burkholderia is an important bacterial genus containing species of ecological, biotechnological, and pathogenic interest. With their taxonomy undergoing constant revision and the phenotypic similarity of several species, correct identification of Burkholderia is difficult. A genetic scheme based on the recA gene has greatly enhanced the identification of Burkholderia cepacia complex species. However, the PCR developed for the latter approach was limited by its specificity for the complex. By alignment of existing and novel Burkholderia recA sequences, we designed new PCR primers and evaluated their specificity by testing a representative panel of Burkholderia strains. PCR followed by restriction fragment length polymorphism analysis of an 869-bp portion of the Burkholderia recA gene was not sufficiently discriminatory. Nucleotide sequencing followed by phylogenetic analysis of this recA fragment differentiated both putative and known Burkholderia species and all members of the B. cepacia complex. In addition, it enabled the design of a Burkholderia genus-specific recA PCR that produced a 385-bp amplicon, the sequence of which was also able to discriminate all species examined. Phylogenetic analysis of 188 novel recA genes enabled clarification of the taxonomic position of several important Burkholderia strains and revealed the presence of four novel B. cepacia complex recA lineages. Although the recA phylogeny could not be used as a means to differentiate B. cepacia complex strains recovered from clinical infection versus the natural environment, it did facilitate the identification of clonal strain types of B. cepacia, B. stabilis, and B. ambifaria capable of residing in both niches.     

Bayesian modeling of recombination events in bacterial populations

Background  

We consider the discovery of recombinant segments jointly with their origins within multilocus DNA sequences from bacteria representing heterogeneous populations of fairly closely related species. The currently available methods for recombination detection capable of probabilistic characterization of uncertainty have a limited applicability in practice as the number of strains in a data set increases.     

Site-directed mutagenesis and virulence assessment of the katG gene of Mycobacterium intracellulare

Mycobacterial catalases have been suggested as acting as virulence factors by protecting intracellular mycobacteria from reactive oxidative metabolites produced by host phagocytes. Mycobacterium intracellulare , like many other mycobacteria, produces two proteins with catalase activity: a heat-stable catalase (KatE) and an inducible, heat-labile catalase peroxidase (KatG). The M. intracellulare katG gene was cloned, and a plasmid derivative with a 4 bp insertion in the katG coding sequence was constructed and used for site-directed mutagenesis of M. intracellulare 1403 (ATCC 35761). The resulting katG mutant was highly resistant to isoniazid (INH), showed an increased sensitivity to H₂O₂ and had lost peroxidase and heat-sensitive catalase activity but retained heat-stable catalase activity. The plasmid carrying the katG frameshift allele was also used for mutagenesis of the mouse virulent M. intracellulare isolate D673. After intravenous injection into BALB/c mice, D673 and the isogenic katG mutant showed the same growth kinetics in the spleen, liver and lungs of the infected mice. Our results demonstrate that the KatG catalase peroxidase mediates resistance to H₂O₂ and susceptibility to INH but is not an essential virulence factor for the survival and growth of M. intracellulare in the mouse.     

Minimum taxonomic criteria for bacterial genome sequence depositions and announcements

Multiple bioinformatic methods are available to analyse the information encoded within the complete genome sequence of a bacterium and accurately assign its species status or nearest phylogenetic neighbour. However, it is clear that even now in what is the third decade of bacterial genomics, taxonomically incorrect genome sequence depositions are still being made. We outline a simple scheme of bioinformatic analysis and a set of minimum criteria that should be applied to all bacterial genomic data to ensure that they are accurately assigned to the species or genus level prior to database deposition. To illustrate the utility of the bioinformatic workflow, we analysed the recently deposited genome sequence of Lactobacillus acidophilus 30SC and demonstrated that this DNA was in fact derived from a strain of Lactobacillus amylovorus. Using these methods researchers can ensure that the taxonomic accuracy of genome sequence depositions is maintained within the ever increasing nucleic acid datasets.